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ABSTRACT 

We present a catalog of 2,990 extended sources in a l°xl° area centered on 

K^ , M33 using the MegaCam camera on the 3.6m Canada-France-Hawaii telescope 

(N \ (CFHT). The catalog includes 599 new candidate stellar clusters, 204 previously 

Q \ confirmed clusters, 1969 likely background galaxies and 218 unknown extended 

objects. We present ugriz integrated magnitudes of the candidates and confirmed 

^ \ star clusters as well as full width at half maximum, ellipticity and stellarity. Based 

O ! oil the properties of the confirmed star clusters, we select a sub-sample of highly 

probable clusters composed of 246 objects. The integrated photometry of the 

complete cluster catalog reveals a wide range of colors from -0.4 < (g-r) < 1.5 

rS \ and -1.0 < (r-i) < 1.0 with no obvious cluster subpopulations. Comparisons 

c3 \ with models of simple stellar populations suggest a large range of ages some as 

old as ~10 Gyrs. In addition, we find a sequence in the color-color diagrams 

that deviates from the expected direction of evolution. This feature could be 

associated with very young clusters (< lO^yrs) possessing significant nebular 

emission. Analysis of the radial density distribution suggests that the cluster 

system of M33 has suffered from significant depletion possibly due to interactions 

with M31. We also detect a gap in the cluster distribution in the color-color 

diagram at (g-r) ~ 0.3 and (u-g) ^ 0.8. This gap could be interpreted as an 

evolutionary effect. This complete catalog provides promising targets for deep 

photometry and high resolution spectroscopy to study the structure and star 

formation history of M33. 



-2- 
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Introduction 



Star clusters provide important information for understanding the formation and evolu- 
tion of galaxies. Such systems are useful for highlighting substructures of their host galaxies 
and for revealing their merging history. In particular, the galaxy formation process can be 
traced through the ages, metallicities, and kinematics of star clusters. Due to their prox- 
imity, galaxies in the Local Group provide us with ideal targets for detailed studies of star 
cluster properties. While the star cluster systems of the Milky Way and M31 have received 
close attention, the third spiral galaxy in the Local Group, M3 3, has been less studied. 
At a distance of 870 kpc (distance modulus = 24.69; iGalleti et al.l 120041 ). M33 is the only 
nearby late-type spira l galaxy (Scd). With a large angular size and inclination of i=56° 
(jRegan &: Vogellll994l ). M33 is a suitable galaxy for studies of its stellar constituents. 



The re have been a number of M33 cluster catalogs published since the pioneering 
work of Hiltnei (|l960 ). An extensive and complete catalog can be found in the work of 
Sarajedini fc Mancond (120071 . hereafter SM), which merged all of the modern catalogs com- 
piled before 2007. The catalogs that have appeared after the publication of SM have been 
incorporated into the web-based version of the SM catalog^. This updated version of the 
catalog contains 595 candidates of which 349 are confirmed clusters based on Hubble Space 
Telescope (HST) a nd high-resolution groun d-based imaging. The most recent work in this 
field corresponds to IZloczewski et al.l ( l2008l . hereafter ZKH) using the MegaCam camera on 
CFHT. This study presents a catalog of 4,780 extended sources in a 1 deg^ region around 
M33 which includes 3,554 new candidate stellar clusters. 



As pointed out by SM, the sample of clusters in M33 suffers from significant incom- 
pleteness. While HST and its several instruments have been successfully used in the search 



""^Based on observations obtained with MegaPrime/MegaCam, a joint project of CFHT and 
CEA/DAPNIA, at the Canada-France-Hawaii Telescope (CFHT) which is operated by the National Re- 
search Council (NRC) of Canada, the Institut National des Science de I'Univers of the Centre National de 
la Recherche Scientifique (CNRS) of France, and the University of Hawaii. 

"■http://www. astro. ufl.edu/^ata/cgi-bin/ni33_cluster_catalog/index.cgi 
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for star clusters, the small field of view permits surveys only over a limited region of the 

2 



galaxy. With the most recent contribution of Z KH, this area has 



)een increased to 1 deg' 



centered on M33. However, as we discussed in ISan Roman et al.l (120091 ). the ZKH catalog 
has largely overestimated the number of clusters, due to a possible systematic misidentifica- 
tion, where only around 40% of the 3554 proposed candidates are likely to be actual stellar 
clusters. For these reasons we have undertaken the present study. This paper is organized 
as follows: Section 2 describes the observations and data reduction while section 3 discusses 
the adopted search method and the integrated photometry of the clusters. The analysis of 
the photometric properties and comparison with other galaxies are in Section 4. Finally, 
Section 5 presents a summary. 



2. Observations and Data Reduction 

The observations for the present study were obtained using the Queue Service Observ- 
ing mode at the 3.6 m Canada France Hawaii Telescope (CFHT). The data are available 
on-line through The Canadian Astronom y Data Centre archiv e and were obtained as part 



of "The M33 CFHT Variability Survey" ( iHartman et al.ll2006[ ). The images were taken us- 



ing the MegaCam/MegaPrime wide-field mosaic imager which contains 36 individual CCDs 
that combine to offer nearly a full l°x 1° field of view with a high angular resolution of 
0.187"pixel~^. MegaCam operates with a set of g'r'i'z' filters very similar to those of the 
Sloan Digital Sky Survey (SDSS) but a slightly different near-UV filter called u*. This filter 
was designed to maximize the capabilities of the instrument at short wavelengths and its 
effective wavelength is ~ 200 A redder than the standard u' filter. All of the archival images 
were pre-processed by the CFHT's Elixir project. This pipeline includes the standard steps 
of overscan and bias subtraction, fiat-fielding, fringe correction, masking of bad pixels and 
merging of amplifiers. The Elixir project also provides a preliminary photometric calibration 
for each image. 

In order to facilitate the search for cluster candidates, only the best available images 
were analyzed consisting of 15 u*, 15 g', 14 r', 28 i' and 3 z'. Median seeing values of all 
analyzed images are < 0.8" in g', r' and i' filters and ~ 0.6" in u* and z'. Prior to the data 
analysis, each field of 36 individual CCDs were combined into a single master image. This 
process was done using the software module Swarp 2.16.14 of the Terapix pipeline which is 
mainly dedicated to the processing of MegaCam data. This specific module involves resam- 
pling of the individual images as well as co-adding the diffe rent exposures in an optimum 



way so that the point-spread function (PSF) is not distorted (JBertin et al.ll2002l . for details) 
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Each final combined master image was divided into two sub-fields, including an overlapping 
area, to deal with the spatial variability of the PSF. 



The M33 images are extremely crowded making the construction of point-spread func- 
tions quite challenging. In order to per form accurate standard profile-fitting photometry, we 
used DAOPHOT/ALLSTAR routines (1Stetsonlll994l ) in an iterative way. First, we found all 
of the stars on each image and produced small-aperture photometry for them. We then used 
the DAOPHOT/PICK routine to select a set of 1000 reasonable candidates to be used as PSF 
stars. After deleting those with bad pixels nearby, we subtracted the stars with surrounding 
neighbors to help isolate the PSF stars. The resulting list of more than 500 stars, in all cases, 
was used to create a PSF for each of the images. The shape of the PSF was made to vary 
quadratically with position on the frame. To improve the PSF, we created an image where all 
the neighbors and stars that do not fit the first PSF were subtracted, obtaining an improved 
second-generation PSF over this subtracted image. Appropriate aperture corrections were 
calculated from isolated unsaturated bright stars with photometric errors smaller than 0.01 
mag. Since the correction varies with radius from the center of the images, a polynomial fit 
was applied to the aperture corrections in order to obtain the final instrumental magnitudes. 
All frames were matched u sing DAOMATCH/DA OMASTER routines to obtain common 
stars in all filters. Following iHartman et al.l (|2006[ ). the photometric calibration provided by 
the Elixir pipeline was applied using the zero-point values. I n addition and to deal with the 
differences between u* and u', we applied the equations from I Clem et al.l (j2008[ ) to transform 
the photometry from u*g'r'i' z' to u'g'r'i' z' . 



The integrated magnitudes and colors for each candi date cluster h ave been calculated 
using the aperture ph otometry routines in DAOPHOT (IStetsonI 119871). To be consistent 



with previous authors ( ISan Roman et al.l l2009t ISarajedini et al. 



2007 



Chandar et al. 1999 



200ll ). we have adopted an aperture radius of 2.2" for the magnitude measurements and 
1.5" for the colors. The background sky is always determined in an annulus with an inner 
radius of 3.5" and an outer radius of 5.0". No aperture corrections have been applied to 
the extended objects, such as the star cluster candidates. Once again, these magnitudes 
have been photometrically calibrated to the SDSS standard system. To derive accurate 
positions of the clusters and to estimate properties such as ellipticity and full width at half 
maximum (FWHM), we have applied the Sextractor v2.5.0 ( jBertin fc Arnoutslll996l ) image 
classification algorithm. 



-5- 

3. A New Catalog of Star Cluster Candidates in M33 

3.1. Cluster Search Method 

Our detection method is based on the fact that at the distance of M33, non-stellar 
objects are expected to be more extended than the PSF. After subtracting the stellar PSF 
from all of the sources in our frames, extended objects leave a doughnut-shaped appearance, 
as they are under-subtra cted in the wi ngs and over-subtracted in the center. We have used 



DAOPHOT/ALLSTAR (IStetsodll994l ) to produce residual images free of all PSF sources. 
We have trained our eyes to recognize the residual pattern of candidate clusters. While most 
background galaxies show either a spiral arm structure or an elongated pattern, the candi- 
date stellar clusters show some level of assembly. As an illustration. Fig. [1] shows original 
and residual images of different types of extended objects. After visual inspection of the 
residual images as well as analysis of the original ones, this technique leaves us with a total 
of 2,990 extended objects: 803 candidate clusters, 1,969 galaxies and 218 unknown objects. 
From the total number of candidate clusters, 204 were previously identified clusters in the 
SM updated website and considered confirmed clusters based on HST and high-resolution 
ground-based imaging. 



The 12 ACS/HST fields examined in ISan Roman et al.l (J2009[ ) included 72 of the can 



didate clusters in the present catalog, where 51 turned out to be genuine star clusters. This 
suggests that around ~ 70% of the proposed candidates will be actual stellar clusters. How- 
ever, from the 349 guaranteed clusters listed in the updated SM catalog, our catalog only 
recovers 204 objects implying missing objects mostly in the center of the galaxy, not surpris- 
ing since the method is less effective in extremely crowded regions. 

Comparison with the similar study of ZKH reveals unexpected discrepancies. From 
the total 1,752 common objects between both catalogs, on l y 124 sources were classified as 



candidate clusters by both authors. As ISan Roman et al.l ( l2009l ) argue, the total number 



of 'true' cluster candidates in the ZKH catalog is not likely to be larger than ~ 40%. This 
suggests a systematic misidentification in the candidate object pattern or a defective PSF 
subtraction. Fig. [2] shows the photometric differences between the two studies. We find a 
mean difference of <Ag> = 0.15 ± 0.02, <Ar> = 0.10 ± 0.02 and <Ai> = 0.04 ± 0.03 
while the offsets for the colors are <A(g-r)> = 0.018 ± 0.012 and <A(r-i)> = 0.037 ± 0.030. 
The disagreement in the magnitude offsets disappears in the color offsets, which indicates 
that the photometric variation corresponds to the different adopted apertures in each study. 
These photometric differences are not unexpected for integrated photometry of extended 
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objects such as star clusters (ISan Roman et al.l 120091 ) 



3.2. Highly Probable Clusters 

We have obtained stellarity, full width at half- maximum (FWHM), and ellipticity for 
the total sample of extended sources by applying the Sextractor software to the target im- 
ages. We have compared the photometric parameters of the known M33 star clusters in 
our sample with the parameters of candidate objects to find a suitable criterion to select 
highly probable clusters. Fig. [3] shows the distribution of the three Sextractor parameters. 
The stellarity parameter object classification of Sextractor allows us to examine our visual 
object classification with a more systematic algorithm. Based on its definition, a stellarity 
of 1 corresponds to a point source (star) and a stellarity of to a resolved object (galaxy). 
Considering the pixel scale of the CCD {0.1S7" pixel~^), a typical seeing of 0.7" and the 
distance of M33 (~ 870 kpc), a mean cluster size of 4pc will appear in our images as a point 
source object of ~ 0.9" implying stellarities around 1. The distribution of elhpticities peaks 
between e=0.05-0.2 with an extended tail reaching 0.5 in both samples. Panel b) shows the 
normalized FWHM assuming a mean seeing of 0.7". The FWHM of our sample has two 
peaks: at FWHM ~ 1.5 that agrees with the confirmed cluster distribution and another 
peak at FWHM ~ 1.1 not associated with the confirmed cluster distribution. Fig. [3] c) 
shows that the distribution of stellarity for the confirmed clusters has a strong peak at ~ 1 
with a weak peak < 0.2. This distribution of stellarities suggests that a significant number 
of confirmed star clusters will be missing in the catalog if our main source of classification 
were the Sextractor classification algorithm. Although the stellarity is a very useful detec- 
tion parameter to distinguish between point sources and non-point sources, the possibility 
of extended or partially resolved clusters in our images, means that the stellarity parameter 
must be used cautiously. 

Based on the properties of the confirmed star clusters, we selected a sample of highly 
probable clusters that satisfy the following criteria: a) Ellipticity < 0.4; b) 1.1 < FWHM(px/ 
seeing) < 2.2; and c) stellarity > 0.6. We show in Fig. [3]d) the correlation between FWHM 
and ellipticity where the filled area corresponds to the specified selection criteria. A min- 
imum condition has been applied to the FWHM in order to avoid stellar contamination 
and small stellar associations. No information on color or magnitude was used for selecting 
the candidates. This sub-sample of highly probable clusters contains 246 objects and has 
been designated as class 2. Analysis of the contamination of this subsample using a similar 
technique as above suggests that ~ 85% of highly probable clusters will be genuine stellar 
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clusters. To illustrate the available data, Table [T] shows an excerpt of the complete extended 
source catalog where the last column corresponds to our proposed classification of the ob- 
jects: l=galaxy, O=unknown extended object, l=candidate star cluster, 2=highly probable 
cluster and 3=confirmed cluster based on the SM catalog. Based on this classification, the 
catalog contains 599 new candidate clusters (353 candidate clusters (class 1) and 246 highly- 
probable clusters (class 2)). Tabled] can be found in its entirety in the electronic edition of 
the Journal. The sample of highly probable clusters as well as the guaranteed clusters in SM 
will be used as targets for future follow-up imaging and spectroscopic observations. 



4. Analysis 
4.1. Color-Magnitude Diagrams 



Intrinsic properties such as age, metallicity, and reddening govern the integrated mag- 
nitudes and colors of clusters. As described in the previous sections, we have performed 
aperture photometry of the candidate (class 1 a nd 2) and con fi rmed (class 3) star clusters. 
In addition, we have made use of the equations in iTucker et al.l ( 120061 ) to transform our pho- 
tometry into the SDSS ugriz standard filters. Fig. H] shows the color-magnitude diagrams 
(CMDs) and color distributions of our sample (class 1 and 2) as compared with the confirmed 
star clusters from the SM catalog (class 3). The magnitude distribution of our sample con- 
tains more faint star clusters than SM. The faintest clusters reach g ^ 22 that corresp onds 



to Mg ~ -3.0, assuming a distance modulus of ( m — M)n= 



=24.69 ( 



Galleti et al.l 120041 ) and 



an average reddening correction of E(V-I)=0.06 ( Sarajedini et al.ll2000l ). The color range 
of our sample is significantly wider than the color range of the confirmed clusters: -0.4 < 
(g-r) < 1.5 and -1.0 < (r-i) < 1.0. The lower panels show a unimodal distribution with a 
strong peak at (g-r) ~ 0.1 and (r-i) ~ 0.2 having extended tails redward in (g-r) color and 
blueward in (r-i). 



4.2. Color-Color Diagrams 



Fig. |5]shows color-color diagrams of the candidate star clusters. To c ompare with simple 



stellar population s (SSP), two d i fferen t sets of models have been used: iBruzual fc Chariot 



(120031, BC03) and iKotuUa et al.l (12009L Galev). BC03 models correspond to an evolutionary 
track for an instantaneous burst and a Salpeter IMF while Galev models correspond to a 
customized set provided to us by Ralf KotuUa and the Galev team. The Galev SSP models 



were run assuming Geneva evolutionary tracks with a minimum age and time resolution of 
0.1 Myrs until 100 Myrs, and a time-step of 1 Myrs for older ages. The models were run with 
a Salpeter IMF (1-120 Mq) for different metallicities. It is important to note that Galev 
models include contributions from nebular emission, considering the continuum nebular and 
also emission lines. All of the cluster s have been shifted by a line-of-sight red dening value of 



E(V-I )=0.06 ( ISarajedini et al.ll2000l ) and adopting an extinction relation from lCardelli et al. 



(J19891 ). 



Comparisons between the integrated cluster colors and the predictions of stellar pop- 
ulation models can provide age estimates that are potentially useful for studies of galaxy 
evolution. However, uncertainties in ages derived from multi-color photometry come not 
only from the photometric errors, but also from reddening corrections and uncertainties in 
the metal abundance of each cluster. Furthermore, the face-on view of the galaxy and the 
numerous spiral arms produce a broad range of reddening in M33 that can scatter the in- 
tegrated colors of individual clusters. Another effect that is important in this regard is the 
dispersion in the integr ated colors due to st ochastic effects and these can vary significantly 



along the age sequence (jGirardi et al.lll995f ). The adopted IMF in the SSP also contributes 
to uncertainties in the models. Given these points, we have not attempted to estimate ages 
based on the integrated photometry of the clusters. In any case, the color-color diagrams 
reveal a number of interesting features. 



A significant fraction of 'bluish' clusters that occupy a unique location in the diagram 
appear in both panels of Fig. [5] at colors (r-i) < -0.2 and (u-i) < 0.8. These clusters repre- 
sent a finger-like feature that deviates from the expected direction of evolution. At least five 
confirmed clusters from SM are associated with the feature, supporting the genuine cluster 
nature of these objects. The age estimates for three of them shows ages ~10^yrs. Based 
on the lower panel in Fig. [5l the 'finger' feature could be associated with the presence of 
a significant population of very young clusters (< lO^yrs) exhibiting nebular emission. The 
position of many of them below the theoretical line is consistent with internal reddening 
shifting their colors to redder values consistent with the dusty clouds in which they are 
born. Fig. O shows the spatial distribution of these very young clusters on a Galex FUV 
image. A close inspection of this image suggests that all the clusters are associated with 
regions of star formation activity. 



The recent work of iGrossi et al.l (J2010l ) analyzes multi-wavelength observations of 32 



young star clusters and associations in M33. The sample was selected from catalogs of emis- 



-9- 



sion line objects based on their round shape and their position in regions of the galaxy that 
are not too crowded in the Ha map. All of the objects have oxygen abundances of 8 < 12 + 
log(0/H) < 8.7 ari d have 24/ini counterparts in the Spitzer/MIPS map. Comparison of the 
Grossi et al.l ( 120101 ) catalog and ours reveals 10 common objects, all of them associated with 
the previously mentioned 'finger' feature. Table [2] provides a cross-identification of the com- 
mon objects and includes the ages, extinctions and reddenings obtained from their spectral 
energy distribution (SED) fitting technique. This result confirms the young age of these clus- 
ters, younger than ~ 12 Myr, and their relatively high extinction, A^,, between 0.5 and 0.9. 
Of the 10 common objects, 2 of them have been previously confirmed as genuine star clusters. 



Prior to this study, the M33 cluster system presented an age range of star clusters be- 
tween 10 Myrs - 10 Gyrs with the majority of clusters with ages around 100 - 400 Myrs. Fig. 
[5] reveals a wider age range of 1 Myrs - 10 Gyrs with at least ~ 50% of objects corresponding 
to young clusters (< 100 Myrs) and ~ 10% of the total corresponding to very young clusters 
with nebular emission (finger-like feature). In spite of the redward tail at (g-r) > 0.7 that 
is probably caused by reddening, the diagram suggests the presence of an old population at 
least as old as 10 Gyr. 

Furthermore, the integrated colors of the nucleus of M33 (RA= 01:33:51.02; DEC= 
30:39:36.68) have been plotted in the color-color diagrams of Fig. \5\ and Fig. O After 
examining the curve of growth for different aperture diameters, we have adopted an aperture 
radius of 4.4" for the magnitude measurements and a background sky of 7.5" and 9" . The 
integrated light shows a blue nucleus with u=14.9, (u-i)=0.92, (u-g)=0.63 and (r-i)=0.33. 



4.3. A Gap in the Diagram 

The color-color diagrams plotted herein reveal a gap in the distribution of s tar clusters 
center ed at (g-r) ~ 0.3 and (u-g) ~ 0.8. A similar anomaly was discovered by iBica et al. 
(Il99ll ) among LMC star clusters in the (U-B) vs. (B-V) diagram. The LMC gap was 
noticed at (U-B) ~ 0.19 and (B-V) ~ 0.47 with an approximate width of 0.1 mag in both 
colors. The upper panel in Fig. [7] shows a small region of the LMC an d M33 color-color 
diagrams for better visualization of the gaps. The LMC data are from iBica et al. ( 19961) 



and h ave been converted to the SDSS ugriz system using relations publish ed by lJester et al. 
( 2005h. Constant reddening values of E(V-I)=0.06 JSaraiedini et al.ll200oh and E(B-V)=0.1 
(JAlves et al.ll2002l ) have been adopted for M33 and LMC, respectively. Although the M33 



gap appears bluer in the diagram, both gaps correspond to a similar range in age. The 
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lower panel in Fig. [7] shows the color distribution of the M33 and LMC clusters. To avoid 
contamination by clusters in the finger-like feature, which corresponds to a different age 
range, only clusters with (u-g) > 0.4 have been considered in the construction of the color 
distribution. The gap in M33 looks smaller and slightly redder than the gap in the LMC. 
The offset could be a consequence of the reddening correction we have applied which does 
not account for dust internal to each galaxy. 



Several authors ( jRenzini fc Buzzonil Il986t ISweigart et al.l Il990l ) have interpreted the 
LMC gap as being produced by the red giant branch phase transition. This transition 
would be pro d uced by stars at the helium flash stage and would flt theoretical predictions. 
Girardi et al.l ( 119951 ) disagrees with this interpretation and argues that the lack of clusters 
in this region is determined by the natural dispersion of the colors. The stochastic effects 
on the mass distribution of stars could produce the dispersion in the colors so no additional 
peculiarities would be needed in the stellar models in order to reproduce this feature of the 
diagram. However, the discovery of the gap among M33 clusters supports the presence of 
an evolutionary effect at that particular age as the origin of both gaps. 



4.4. Additional Comparisons With Other Galaxies 

Fig. [H] presents the color-magnitude diagrams of the cluster system of our sample, 
the Mil ky Way (MW), M 31 and the Large MagellanicCloud (LMC). The MW data be- 
long to iLata et al.l ( 12002| ) (open clusters) and JHarrisi ( 1l996l ) (globular cluste r s). T he M31 
data correspond to the candidate and confir r ned cl usters in jPeacock et al.l ( l2009l ). Data 
for the LMC cluster system from JBica et al.l ( 1l996l ) have also been plotted. In order to 
compare the different cluster systems, we have plotted absolute magnitudes and redden- 
ing corrected colors. The MW absolute magnitudes are taken directly from the above- 
mentioned catalogs assuming a speciflc distance m odulus and reddening for each cluster. We 
have adopted a LMC distance modulus of 18.50 (JFreedman et al.ll200ll) and 24.36 for M3 1 



f lVilardellet"aDl2010h. C onstant reddening values of E(V-I)=0.06 JSaraiedini et al.ll200ol ) 



E(V-I)=0.1 fJDurrell et al.ll200ll ) and E(B-V)=0.1 fJAlves et al.ll2002h have been adopted for 
M33, M31 and the LMC, respectively. If needed, we have used the Jester et al (2005) trans- 
formations to convert absolute magnitudes into the g-band. The dashed lines represent the 
division of Galactic globular clusters at (B-V)o = 0.5. 



No distinct cluster subpopulations can be identifled within the M33 cluster system like in 
the LMC or MW. However, the integrated colors of the very young clusters are not necessar- 
ily a reflection of their ages because they could be affected by nebular emission. Considering 
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the significant number of tliis type of object in our sample, tfie color distribution could be 
distorted and appear unimodal when in fact it is not. When the nebular emission clusters 
are removed from the analysis (see Fig. [7]), the color-magnitude diagram shows a possible 
bimodality. The M33 and LMC systems are dominated by blue clusters, (B-V)o < 0.5, in 
contrast with the redder M31 system. However, while the red cluster subpopulation of LMC 
occupies a very narrow (g-r)o region, M33 red clusters populate a significantly wider color 
range more similar to M31 red clusters. 



When comparing the absolute magnitudes of the cluster systems, we see that the bright- 
est clusters in the MW and LMC reach luminosities of M^ ~ -9.5; however, the brightest 
clusters in M33 correspond to M^ ~ -8, more than one magnitude fainter. This effect could 
be explained by the relation between the star formation ra te (SFR) of a galaxy and the 
maximum mass/luminosity of its star clusters ( jLarsenl 120021 ). The empirical relation sug- 
gests that galaxies with high SFRs form proportionally more clusters, and as a consequence, 
the cluster mass function reaches higher masses. Assuming a SFR of ~ 0.45 M iy, yrs "-^ 
the cl uster system of M33 would fit onto this relation reasonably well (see Fig. 1 



feooshl. With the slightly h igher SFR for the MW, LMC and M31 (eg. IXang et al. 



Bastian 



(120091), 



Robitaille &: Whitneyl ( 120101 )) these systems will produce brighter clusters. In addition, en- 
vironmental variations, such as the mass/lur ninosity of the gala xy, can play a role in the 
color- magnitude diagram of a cluster system (jMieske et al.ll2010[ l. 



Fig. M presents the color-color diagrams of the M33 cluster system using our sample, 
M31 and LMC. The sources of the data are the same as those given above. As a reference, 
SSP models from the Galev team (JKotuUa et al.l 120091 ) with a metallicity of z=0.0004 have 
been overplotted. To identify different time periods, the star symbols correspond to 10^, 10^, 
10^, 10^ and 10^° yrs. The same constant reddening value has been adopted for each sample 
as in the previous figure. 



The wide color range of the M33 clusters, -0.4 < (g-r) < 1.5, overlaps entirely with 
the young- intermediate age system of the LMC and with the older M31 system. The broad 
range of colors implies a large range of ages, suggesting a prolonged epoch of formation. 
Based on this evidence, the majority of the clusters will be young-intermediate age objects 
although we would expect clusters older than 10 Gyrs. The diagram also shows that a small 
group of M31 clusters occupies the unique area of the 'finger' feature, however the region 
seems to be significantly more populated in M33 than in these two galaxies. 
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When comparing M33 with similar morphological type galaxies such as NGC 300 or 
MlOl, M33 seems to posses a unique very young star cluster population. The color dis- 
tribution of candidate clusters in MlOl is similar to M33 candidate clusters but no evi- 



dence of very young clusters with nebular emission has been found (JBarmby et al.l 12006 



Chandar et al.ll2004j ). Although NGC 300 is nearly a twin galaxy of M33 in terra s of Hubble 
type and mass, there are several differences between them (JGogarten et al.ll2010l). NGC 300 
appears to have globular clusters similar to those of the Milky Way (INantais et al.l 120101) 



and a metallicity gradient consistent with stars formed prior to 6 Gyrs ago ( jGogarten et al. 



20101 ). Environmental factors may play a key role in the star formation history of M33, 
as NGC 300 is isolated from other galaxies w hi le M33 appears to be inte racting with M31 
( JMcConnachie et al.ll2009l : iPutman et al.ll2009l ). iMcConnachie et al.l (|2009[ ) propose a plausi- 
ble M31-M33 interaction model that reproduces with good agreement the observed distances, 
angular positions and radial velocities of these galaxies as well as the well-known HI warp 
in M33. In this simulation, M33 starts its orbit around M31 about 3.4 Gyrs ago reaching 
pericenter (r ~ 56 kpc) around 2.6 Gyrs ago. After it passes apocenter (r ~ 264 kpc) about 
900 Myrs ago, M33 would be approaching M31. This close encounter could have triggered 
an epoch of star formation in M33. The significant population of very young clusters with 
nebular emission and their association with star formation regions are evidence in support 
of recent star formation activity in M33. 



Many stu dies have shown that interacting/merger environments form large populations 



of clusters (e.g JWhitmore fc Schweizerl ( 119951 ): iGallagher et al.l (l200ll )). especially very young 



clusters. We would expect to see very young clusters still embedded in their dust cocoons in 
these disturbed systems. Yet, their color-color diagrams do not exhibit as prominent a finger- 
like featur e due to nebular emission aro und very young clus ters as compared with M33 (see 



Antenna ( IWhitmore fc Schweizeiill995l ): Stephan's quintet (IGallagher et al. 



2001)). In the 



unusual environment of Hickson compact group 31 (HGC 31). IGallagher et al.l (120101 ) found 
a large population of < 10 Myr star clusters with strong nebular emission, similar to the one 
found in the present study. The main galaxies that make up HGC 31 are disrupted under 
the presence of strong gravitational interactions and show tidal structures. The star cluster 
candidates with nebular emission appear throughout HCG31, specifically concentrated in 
the interaction regions. The existence of these very young star clusters seems to be the 
consequence of active recent and ongoing star formation in HGC 31. The similarities between 
the M33 cluster system and that of HGC 31, which is a strongly interacting environment, 
support two important assertions. First, the finger-like feature is a genuine characteristic 
and not an artificial effect due to the contamination of our cluster sample. Second, the past 
interactions between M33 and M31 have likely had significant impact on the properties of 
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the M33 cluster system, especially the youngest clusters. 



4.5. Spatial Distribution 

In order to analyze the spatial distribution of clusters with different ages, we have di- 
vided our sample into two groups based on comparisons with SSP models. Based on BC03 
models, clusters with (r-i)o ~ and (g-r)o ~ 0.1 have ages of ~ lO^yrs. We are going to 
consider clusters with (r-i)o > and (g-r)o > 0.1 as red or old clusters. The remaining clus- 
ters we categorize as blue or young objects. This partition minimizes the contamination of 
the old (red) clusters by the young clusters in the finger-like feature which exhibit integrated 
colors that are redder than expected. In order to place the cluster density d istribution in 



the context of the field stars, we have made use of the iHartman et al.l (120061 ) star catalog 



constructed over the same MegaCam/CFHT images used in the present study. 

Fig. [To] shows the cumulative radial distributions of the young and old cluster popula- 
tions as compared with the blue (young) and red (old) field star populations. Blue clusters 
follow a spatial distribution similar to the blue field stars. The distribution also suggests 
that younger (bluer) clusters are more centrally concentrated as compared with older (red- 
der) clusters. The red clusters are more dispersed in a wider region than the bluer ones, 
indicating that the majority of the red (old) clusters likely belong to the halo while the bluer 
(younger) clusters generally belong to the disk of M33. Analysis of these distributions using 
a Kolmogorov-Smirnov (K-S) test shows that there is a greater than 99.9% chance that the 
old cluster population is significantly different than the young cluster population. 

Fig. [11] shows the radial density distribution of our entire cluster sample. The filled 
circles show the cluster density profile versus deprojected radius, assuming our adopted dis- 
tance modulus of (m — M)o=24.69, while the open circles show the confirmed clusters from 
SM for comparison. The small dots correspond with the radial density distribution of the 
field stars where the solid line represents the best polynomial fit. The star density distribu- 
tion has been scaled to match the cluster density in the region between R=0.6 - 2 kpc where 
both distributions are likely to have similar completeness levels. 

Inside ~ 0.8 kpc, the cluster profile presents a decrease in density, suggesting some level 
of incompleteness. The cluster profile outside of ~ 2.5 kpc could be reproduced by a power- 
law where the most distant clusters are located at ~ 15 kpc (63 arcmin) from the center of 
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the galaxy. Confirmed clusters from SM show that the M33 cluster system seems to be more 
centrally concentra ted than the field stars, however no other galaxy has been found with 



this characteristic ( jForbes et al.lll996l ). Our radial profiles, shown in Fig. [TTl have reduced 
the discrepancy between the clusters and field stars but the former are still more centrally 
concentrated than the latter as shown by the cumulative distributions in Fig. [TOl The 
pronounced decrease at ~ 15 kpc in the cluster and field-star density distributions suggest 
that this distance may represent the outer edge of both distributions. For a given radial bin in 
the outer region of the galaxy, the density of clusters is significantly lower than the density 
of stars. We note the possibility that the cluster and stellar samples may have different 
completeness properties. In order to minimize the potential impact of incompleteness, we 
restrict the comparison of these samples to the region outside ~ 0.8 kpc from the center of 
the galaxy. If the incompleteness of our sample is the reason for the differences between the 
cluster and field star radial profiles, then we would expect a random bias or perhaps a larger 
incompleteness toward the center of the galaxy. However, the analysis shows a significantly 
lower density of clusters between R = 3 - 9 kpc than in the inner region between R=0.8 - 
3 kpc. 

The ratio of stars to clusters is determined not only by the formation processes but 
also by the destruction processes. If we assume that the formation of star clusters and 
the formation of stars in a galaxy are correlated, then Fig. [11] suggests that the cluster 
system in M33 has suffered from destruction or depletion of clusters at specific radii. Tidal 
interactions when passing through the disk or near massive objects such as giant molecular 
clouds could produ c e tidal shocks that l ead to the ultimate destruction of a cluster (e.g. 



Gieles et al.l (120061 ): iLamers et al.l ( 120051 )). Analysis of the dynamical evolution of these 
clusters is needed to reveal the level of influence of these interactions in the disruption 
process. Other environment effects such as interactions between M33 and M31 can also 
play a role in the depletion or disruption of clusters at preferred galactocentric distances. 



Huxor et al.l (120091 ) discovered the presence of four new outlying star clusters in M33 which 
have large projected radii of 38 - 113 arcmin (9.6 - 28.5 kpc). Based on the asymmetry in 
the distribution of these outer clusters, they suggest the possibility that interactions with 
M31 may have dramatically affected the population of M33 star clusters. Regardless of the 
source of this anomaly, we would need an additional ~ 350 clusters between R=3 - 9 kpc in 
order to match the stellar density in the same region of the galaxy. 

If we rescale the density of stars to match the cluster density in the outer region, 
a notable excess of clusters occurs at R < 4 kpc. This scenario is highly unlikely since 
dynamical destruction processes are more effective near the central region of a galaxy. The 
short lifetime of such a young sample of clusters also makes the cluster migration scenario 
implausible. No case has been found in which the cluster density exceeds the star density 
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in the inner region of a galaxy. Future follow-ups of this sample will test the validity of the 
depletion phenomenon that could have widespread repercussions for our understanding of 
M33's formation and evolution. 



5. Summary 

We present a wide-field photometric survey of M33 extended objects using CFHT/ 
MegaCam images. The resultant catalog contains 2,990 extended sources, including 599 new 
candidate stellar clusters and 204 previously identified clusters. We have investigated the 
photometric properties of the cluster sample, performing ugriz integrated photometry and 
using their morphological parameters. Based on the properties of confirmed star clusters, we 
select a sub-sample of 246 highly probable objects. Analysis of multicolor photometry of the 
candidate clusters reveals a wide range of colors including a finger-like feature in the color- 
color diagrams that deviates from the expected direction of evolution. Color distributions 
of the cluster sample reveal a unimodal distribution. A comparison of the radial density 
distribution for the field stars and our cluster sample suggests that the M33 cluster system 
suffers from a depletion of clusters at all radii. Color-color diagrams also reveal a gap in the 
distribution of star clusters similar to the gap detected among LMC clusters. 

We thank Ralf Kotulla for generously provide us with the customized Galev models 
used in this work. We also appreciate the useful comments of Rupali Chandar. We are 
grateful for support from the United States National Science Foundation via grant number 
AST-0707277. 
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Table 1. Extended Source Catalog 



Id R.A. (J2000.0) Decl. (J2000.0) 



(u-g) (g-r) (g-i) (g-z) EUipticity FWHM (") Stellarity Alter. Id'' Classification'' 



1 


1 31 33.13 


31 04 05.65 


2 


1 31 33.39 


30 38 07.97 


3 


1 31 33.88 


30 58 40.34 


4 


1 31 33.96 


31 09 09.60 


5 


1 31 34.05 


31 07 12.55 


6 


1 31 34.07 


30 40 58.22 


7 


1 31 34.21 


30 36 09.61 


8 


1 31 34.36 


31 00 37.44 


9 


1 31 34.52 


30 53 13.78 


10 


1 31 34.62 


30 39 10.50 


11 


1 31 34.67 


30 20 06.59 


12 


1 31 34.78 


31 01 35.22 


13 


1 31 35.13 


30 17 42.01 


14 


1 31 35.20 


30 48 11.67 


15 


1 31 35.28 


30 32 39.49 


16 


1 31 35.49 


30 45 21.35 


17 


1 31 35.57 


30 27 44.69 


18 


1 31 35.59 


30 27 55.07 


19 


1 31 35.67 


31 00 10.35 


20 


1 31 35.69 


31 07 25.79 


21 


1 31 35.70 


30 17 19.51 


22 


1 31 35.71 


30 36 09.99 


23 


1 31 35.89 


30 52 54.88 


24 


1 31 35.98 


30 23 53.63 


25 


1 31 36.03 


31 05 26.83 


26 


1 31 36.05 


30 53 42.02 


27 


1 31 36.16 


31 06 04.95 



18 



.464 



20 



107 



950 



957 



750 



453 



932 



625 0, 



704 



0.23 


0.88 


0.98 


0.53 


1.83 


0.02 


0.37 


1.58 


0.01 


0.19 


1.46 


0.22 


0.07 


1.21 


0.07 


0.06 


0.83 


0.93 


0.42 


1.04 


0.42 


0.13 


1.11 


0.53 


0.39 


0.00 


0.35 


0.38 


1.57 


0.00 


0.16 


0.86 


0.98 


0.17 


1.01 


0.67 


0.77 


1.07 


0.05 


0.28 


0.99 


0.09 


0.45 


1.28 


0.01 


0.21 


1.33 


0.00 


0.24 


1.11 


0.56 


0.31 


1.44 


0.04 


0.24 


1.28 


0.11 


0.26 


0.95 


0.04 


0.19 


0.89 


0.93 


0.15 


1.24 


0.07 


0.22 


1.22 


0.38 


0.28 


2.57 


0.06 


0.43 


1.55 


0.00 


0.32 


1.29 


0.29 


0.21 


0.89 


0.89 



-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 

2 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
2 
-1 
-1 
-1 
-1 
-1 
-1 



CO 



Note. — Table 1 is published in its entirety in the electronic edition of the Journal. The complete table also includes the original CFHT filters (u*g'r'i'z') 
magnitudes and errors for the candidate star clusters (SC). 
Units of RA are hours, minutes, and seconds, and units of Dec are degrees, arcminutes, and arcseconds. 

^Identification Number in lSaraiedini fc Mancong 1(2003) 

'^Proposed classification: -1 galaxy, unknown extended object, 1 candidate SC, 2 highly probable SC, and 3 confirmed SC (based on SM catalog) 
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Table 2. Cross Identification with Grossi et al. (2010) 



Id (Us) 


Id (Grossi) 


log Age'^ 


A,^ 


E(B-V)^ 


Glassification'^ 


1742 


VGHC 2-84 


6.51±0.14 


0.58±0.11 


0.26±0.03 


3 


735 


C400 


6.65±0.15 


0.52±0.17 


0.23±0.07 


1 


1144 


C129a 


6.35±0.21 


0.93±0.09 


0.39±0.02 


1 


1084 


C121 


6.36±0.13 


0.93±0.06 


0.41±0.01 


2 


970 


B0221 


6.99±0.04 


0.40±0.02 


0.16±0.01 


2 


760 


LOG HII 3 


6.43±0.17 


0.54±0.08 


0.22±0.02 


2 


847 


B0261 


6.53±0.16 


0.56±0.10 


0.24±0.04 


2 


983 


MGM00Em24 


6.37±0.09 


0.61±0.03 


0.28±0.01 


2 


1586 


B0013C 


7.16±0.11 


0.56±0.11 


0.27±0.09 


3 


952 


G403 


6.49±0.14 


0.88±0.11 


0.38±0.03 


2 



'^Gluster properties from Grossi et al. 2010. 
Proposed classification as in Table 1. 
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Fig. 1. — Illustration of the search method adopted. First column corresponds to the original 
images while the second column corresponds to the residual images after PSF subtraction. 
Rows from top to bottom show: a background galaxy, a confirmed star cluster, and two new 
candidate clusters. Each image is shown with the same gray-scale intensity and 15" on a 
side, with north up and east to the left. 



-22- 






I 
"3* 



I 



1.0 
0.5 
0.0 

1.0 
0.5 
0.0 

-0.5 

-1.0 
1.0 

0.5 

0.0 

-0.5 

-1.0 

-1.5 
0.4 

0.2 
0.0 

-0.2 

-0.4 

1.5 

1.0 

0.5 

0.0 

-0.5 

-1.0 

-1.5 



• 
• 
• 


• 
• 
• 

•J •'•• . • 
ijt-\iX^> -i^ ....^.. .'..... 


■ ■ •; .«'*T(t ■ ■ ■•■»;•»/■ .■■•■■ ■ i 


• 

1 1* 1 1 


• 

1 • 1 1 


• • 

• • 
• 

• 


•• 

• • 

• • 

• 


• • ' •* 

• . .•.•••*. *• • • 

-••■» v:.^-^?.-t.-^:i-v:--- ■ 

• : • • •; •:• • 

. •• • . . : 


'■■.■■■■3ft--'.vV •■■■^■".■* 

/ . •• • 


• 

• 


• 
• 



18 



19 



20 



21 



-0.5 



0.0 



0.5 



1.0 



1.5 



Fig. 2. — Comparison of the integrated cluste r photometry from the present study with the 
common objects from IZloczewski et aLl ( 120081 ) . 
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Fig. 3. — Sextractor parameter distributions of the 599 new candidate clusters (open circles) 
and the 204 previously confirmed clusters in the catalog of SM (filled circles). The open 
histograms correspond to the distribution of confirmed clusters while the filled histograms 
correspond to our sample of candidate clusters. The filled area in panel d) shows the selected 
sample of highly probable clusters. 
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Fig. 4. — Color-magnitude diagrams and color distributions of the new star cluster candi- 
dates (open circles/unfilled histogram) as compared with confirmed clusters from SM (filled 
circles/filled histogram) . 
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Fig. 5. — Color-color diagram of the candidate star clusters from the present study (black 
circles) as compared with the confir med clusters from SM (gre en circles). The solid lines 



corre spond with the SSP models of (jBruzual fc CharlotI l2003l . BC03) and (jKotulla et al. 



20091 . Galev) for z=0.019 and z=0. 0004. All the clusters h ave been shifted by a line-of-sight 
reddening value of E(V-I) = 0.06 (ISarajedini et al.l 120001 1 . The red cross corresponds with 
the integrated colors of the nucleus of M33 and the arrows in the top right corners represent 
the direction of the reddening vector. [See the electronic edition of the Journal for a color 
version of this figure.] 
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Fig. 6. — Galex FUV image of M33. The image traces the star formation regions of the galaxy 
through its young stars. The black triangles correspond to the 56 star cluster candidates 
associated with the finger-like feature. [See the electronic edition of the Journal for a color 
version of this figure.] 
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Fig. 7. — Upper panel: Color-color diagram of our sample (filled circles) and LMC (open cir- 
cles) ( JBica et al.lll996l ) in the relevant color ranges. To avoid contamination by clusters that 
exhibit nebular emission, only clusters with (u-g) > 0.4 (dashed line) have been considered 
in the color distribution. Lower panel: Color distribution of our sample (unfilled histogram) 
and LMC (filled histogram). M33 gap can be detected at (g-r) ^ 0.3 and (u-g) ~ 0.8. , and 
the LMC gap at (g-r) ~ 0.3 and (u-g) ~ 1.3. 
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Fig. 8. — Color- magnitude diagrams of star clusters in different galax ies of the Local Group: 
a) The Milky Way: open clusters ( open circles) (Lata et al.l |2002|) and globular clusters 



(filled circles) flHarrislll996f ): b) M31 f lPeacock et al. 



2OO9I ) c) LMC teica et al-lEoOGl ) and c 



M33 confirmed and candidate star clusters from this study. The dashed lines represent the 
division of Galactic globular clusters at (B-V)o = 0.5. 
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Fig. 9. — Color-color diagrams of the candidate star clusters from the present study (black 
circles) as con i pared with the candidate and confirmed clusters in M 31 (green circles ) 
(jPeacock et al.l |2009[ ) and confirmed clusters in the LMC (blue circ les) (JBica et al.l Il996l ) . 
The solid line corresponds to the SSP models of iKotuUa et al.l ( 120091 ) with a metallicity of 
z=0.0004. The star symbols correspond to ages of 10^, 10^, 10^, 10^ and 10^° yrs and the red 
cross corresponds to the integrated colors of the nucleus of M33. \See the electronic edition 
of the Journal for a color version of this figure.] 
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Fig. 10. — Cumulativ e radial distribution for the star clusters (this study) and field stars 
(JHartman et al.ll2006[ ) in M33. \See the electronic edition of the Journal for a color version 
of this figure.] 
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Fig. 11. — Radial density distri bution of our clu s ter sa mple (filled circles) as compared with 
the field stars (small dots) from iHartman et al.l (120061 ) versus deprojected radius. The solid 
line represents the best polynomial fit of the field star radial density. As a comparison with 
previous catalogs, open circles correspond to the confirmed clusters in SM. 



